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Stationary droplets simulated by multi-phase lattice Boltzmann methods lead to spuri- 
ous velocities around them. In this article I report the origin of these spurious velocities 
for one example and show how they can be avoided. 



1. Introduction 

Imagine a stationary droplet in quiescent conditions. For a real system there are, of 
course, no velocities. Yet when this simple system is simulated with a lattice Boltz- 
mann method you will find a flow-field around the drop. Even if you start with a no 
flow initial condition this velocity field will develop. An example of these velocities 
in shown in Figure [l] (a). These velocities are know as "spurious velocities". 

The magnitude of the velocities depends on the details of the method, the ra- 
dius of the drop, the surface tension and the viscosity. There have been several 
studiesEl'BiJ of different methods that tell us about the dependence of the spurious 
velocities on these parameters. In the past the focus of the work has been on re- 
ducing the magnitude of these velocities but to the best of my knowledge no good 
understanding of their origin has been reached. In this article I will explain why we 
see these spurious velocities at all and how, and at what cost, they can be avoided. 

I will first introduce a simple lattice Boltzmann method for model B dynamics 
and show the perfect approach to equilibrium which is free of any spurious currents. 
Then I extend this model to include hydrodynamics and we will see that spurious 
currents suddenly appear. A simple examination shows why there should not be 
any spurious currents according to the continuum equations and I point out the 
terms through which discretization errors drive the spurious currents. I introduced a 
thermodynamically consistent discretization of these terms which leads to a lattice 
Boltzmann method that is free of spurious currents. 
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2. Lattice Boltzmann for Model B 

Model B is a model describing the behavior of binary alloys. We will consider a 
mixture of, say, A and B atoms. The system is characterized by the order parameter 
(f) which represents the difference in the densities of the two components, i.e. (j) — — 1 
represents pure A and (j) = 1 represents pure B. The dynamics of this conserved 
order parameter (j) is then given by 

5F 

dtcji = Vi^V^; ^i=J^ (1) 

where D is a diffusion coefficient and /i is the chemical potential. The system is 
described by the free energy F. For simplicity we will consider a Landau Free 
energy expansion around the critical point which is given by 

F = j + + (2) 
To simulate this equation we use a BGK lattice Boltzmann method given by 

g,(x + V,, i + 1) = .g,(x, t) + -[.g,l'(x, t) - g,(x, t)] (3) 

T 

where the relaxation time r may depend on (j) and the set of discrete velocities 
{vi} will be determined later. The behavior of the model will be determined by the 
moments of the equilibrium distribution g^. We name these moments 

'^ = E5° = E5- F^^gS; Af^J^sSv,. (4) 

i i i i 

A Taylor expansion to second order in the derivatives gives 

(dt + v,V)gO - (dt + v,V)(r ~ + ^,V)g^ = ^[5° - 5.] (5) 

Z T 

We now take the first moment of this expansion and obtain as the equation 

of motion for the order parameter 

dt(j) + VF = V(t - 0.5)dtF + V(t - 0.5) VM (6) 

We have two independent ways to ensure that this equation be equivalent to (|l|). 
Either we choose F = 0, M = /il, t = D + 0.5 or alternatively F = DV/x and 
M — 0. Note that the V9tF term is now third order. In two dimensions we can 
implement both of these approaches with a small velocity set of five velocities 
{vj = {(0, 0), (1,0), (-1, 0), (0, 1), (0, -1)}. 

Both approaches behave equivalently, as expected, but they differ in the ranges 
of stability as will be discussed elsewhereo. Both implementations of this simple 
diffusive lattice Boltzmann algorithm are free from spurious currents. We conclude 
that hydrodynamics is required to see spurious currents. 
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3. Model H binary fluids 

In order to simulate binary fluids we have to introduce a fluid velocity which will 
convect the order parameter (/). This fluid velocity will obey the Navier-Stokes equa- 
tion which is coupled to the order parameter through a thermodynamics pressure 
tensor P'^. We will denote the fluid density as p. The equations of motion area 

pdt(j) + /9V(0u) = VDp{l - 02)v^ + \/D{p - (j)T)yp (7) 
dtp + V(pu) = (8) 
pdtu. + puVu Vjj + VP" + ?7V[Vu + (Vu)^ - tr(Vu)/d] (9) 

where the pressure is p ^ pT ^ r\ is the viscosity and d the number of spatial dimen- 
sions. The pressure tensor P" is defined as 



P" = [(t>d^F -F - K{(f>y^(j> + 0.5V (t)V (/))]! + KV(f>V(f> 



(10) 



We can implement these continuum equations using lattice Boltzmann. With (|^) we 
now define 5j = 4>P- By imposing J^i 9i'^i = 4>P'"-+F and 9^"^i"^i = M +0puu 
we obtain the required order parameter equation with either of the previous 
choices for F and M . The Navier-Stokes equations are obtained with a usual one- 
component lattice Boltzmann method which is coupled to the order parameter 
equation through a body force of VP'^. 





(a) 

Fig. 1. (a) Spurious currents for a drop where the longest velocity vector corresponds to a 
lattice velocity of 3.4 10~^.The density of the matrix is </> = 0.931 and the density of the drop 
is <j> = —1.013. (b) The chemical potential and the spurious diffusion currents for the same drop. 
The longest velocity vector corresponds to a lattice velocity of 2.4 10"'*. The values of ^ vary 
between -1.298 10"^ and -1.181 10"^. 



Doing this leads to the well known spurious velocities as shown in Figure y. We 
see that there are not only spurious currents but also the chemical potential ji is 
no longer constant. Since the implementation of a simple model B did not lead to 
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the spurious velocities it is reasonable to assume that the implementation of the 
additional lattice Boltzmann equation for the total density and the momentum is 
the reason for the occurrence of the spurious velocities. Why should the continuous 
equations in equilibrium be consistent with a quiescent drop without spurious ve- 
locities? This is because the two driving terms VDV/x and VP'^ are related through 

VP" (11) 

This means that a constant chemical potential will lead to zero driving force in 
both equations. 

When we now use this knowledge and replace the driving force VP'^ with (pWn 
we find that the spurious velocities vanish to machine precision. The density of the 
order-parameter in the surrounding fluid is = 0.933 and the density in the drop is 
(j) = -1.015 (cf Figure p. The size of the maximum velocity is now 5.6 10 The 
chemical potential has the value of /i = —1.2103 10~^ and variations are smaller 
than 2 10~^^. So why do the spurious velocities appear in the first place? This 
is because the discretizations of WP^ and 0V/i are very different. We conclude 
that the different discretizations of the driving forces for the order parameter and 
momentum equations are the origin of spurious velocities for lattice Boltzmann. 

There is, however, a caveat. The term 0V/i in its discrete form is not a divergence 
of a scalar field which means that momentum is now only approximately conserved. 
To be able to see the absence of spurious velocities I included a tiny correction term 
in the definition of the momentum that ensures that the total momentum of the 
system does not change. Also my implementation of the algorithm turned out to be 
unstable so I added a small amount of numerical viscosity by multiplying a velocity 
by (1 — 0.001) and adding 0.00025 times the velocity of the four nearest neighbors. 
This rendered the simulations stable. 

4. Conclusions 

I proved that spurious velocities in one particular lattice Boltzmann implemen- 
tation are caused by non-compatible discretizations of the driving forces for the 
order-parameter and momentum equations. The different discretization errors for 
the Forces compete and drive the spurious currents. I believe that the same ar- 
gument holds for all lattice Boltzmann methods that exhibit spurious velocities. 
These spurious velocities can be avoided by ensuring that the discretizations of the 
driving forces are compatible. 
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